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ABSTRACT 

For many years, astronomers have promised that the study of pulsating white dwarfs would ultimately lead 
to useful information about the physics of matter under extreme conditions of temperature and pressure. In this 
paper we finally make good on that promise. Using observational data from the Whole Earth Telescope and a new 
analysis method employing a genetic algorithm, we empirically determine that the central oxygen abundance in 
the helium-atmosphere variable white dwarf GD 358 is 84 ± 3 percent. We use this value to place preliminary 
constraints on the 12 C(a,7) 16 nuclear reaction cross-section. More precise constraints will be possible with 
additional detailed simulations. We also show that the pulsation modes of our best-fit model probe down to 
the inner few percent of the stellar mass. We demonstrate the feasibility of reconstructing the internal chemical 
profiles of white dwarfs from asteroseismological data, and find an oxygen profile for GD 358 that is qualitatively 
similar to recent theoretical calculations. 

Subject headings: methods: numerical — nuclear reactions, nucleosynthesis, abundances — stars:individual 
(GD 358) — stars: interiors — stars: oscillations — white dwarfs 



> 

in 
o 

(N 
O 



^ ■ 

Oh; 

6 ■ 
a: 



The value of the l2 C(a, 7) 16 cross-section at stellar energies 
is presently the result of an extrapolation across eight orders of 
magnitude from laboratory data (Fowler 1986). During the final 
stages of a red giant star, this reaction competes with the triple- 
a process for the available a-particles. As a consequence, the 
final ratio of carbon to oxygen in the core of a white dwarf is 
a measure of the relative rates of the two reactions (Buchmann 
1996). The triple-a reaction is relatively well determined at the 
relevant energies, so the constraint that comes from this mea- 
surement is much more precise than other methods. 

2. FORWARD MODELING 

The initial 3-parameter fits to GD 358 (Metcalfe Nather & 
Winget 2000, hereafter MNW) made it clear that both the cen- 
tral oxygen abundance and the shape of the internal chemical 
profile should be treated as free parameters. We modified our 
code to allow any central oxygen mass fraction (Xo) between 
0.00 and 1.00 with a resolution of 1 percent. To explore dif- 
ferent chemical profiles we fixed Xo to its central value out to 
a fractional mass parameter (q) which varied between 0. 10 and 
0.85 with a resolution of 0.75 percent. From this point, we 
forced Xo to decrease linearly in mass to zero oxygen at the 95 
percent mass point. 

This parameterization is a generalized form of the "steep" 
and "shallow" profiles used in Metcalfe Nather & Winget 
(2000). We originally used these profiles so that our results 
could be easily compared to earlier work by Wood (1990) and 
Bradley Winget & Wood (1993). The latter authors define both 
profiles in their Figure 1. The "shallow" profile corresponds 
approximately to q = 0.5, and "steep" corresponds roughly to 
q = 0.8. However, in our generalized parameterization we have 
moved the point where the oxygen abundance goes to zero from 
a fractional mass of 0.9 out to a fractional mass of 0.95. This 
coincides with the boundary in our models between the self- 
consistent core and the envelope, where we describe the He/C 
transition using a diffusion equilibrium profile from the method 
of Arcoragi & Fontaine (1980) with diffusion exponents of ±3. 
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1. INTRODUCTION 

White dwarf stars are the end-points of stellar evolution for 
the majority of all stars, and their composition and structure can 
tell us about their prior history. Compared to main-sequence 
stars, white dwarfs are relatively simple objects to model. Once 
they pass through the hot planetary nebula nucleus phase, their 
evolution is a slow cooling process largely uncomplicated by 
nuclear burning or other poorly understood effects. Nature has 
been kind by providing several varieties of otherwise normal 
white dwarf stars that undergo continuous and remarkably sta- 
ble pulsations. We can determine the internal structure of these 
pulsating white dwarfs by observing their variations in bright- 
ness over time and then matching these observations with a 
computer model which behaves the same way. 

In the past decade, the observational aspects of white dwarf 
seismology have been addressed by the development of the 
Whole Earth Telescope (Nather et al. 1990). This instrument 
is now mature, and has provided a wealth of seismological 
data on the different varieties of pulsating white dwarf stars. 
In an effort to bring the analytical treatment of these data to 
a comparable level of sophistication, we have recently devel- 
oped an objective global optimization method for white dwarf 
pulsation models based on a genetic algorithm [see §2 of Met- 
calfe Nather & Winget (2000) and references therein for a de- 
tailed discussion of genetic algorithms]. Metcalfe Nather & 
Winget (2000) applied this new approach to observations of the 
helium-atmosphere variable white dwarf GD 358, allowing the 
stellar mass (M»), the effective temperature {T e s), and the mass 
of the surface helium layer (Me e ) to be free parameters. They 
repeated the procedure for several combinations of core compo- 
sition and internal chemical profiles. Among other things, this 
initial study demonstrated that the fits were significantly im- 
proved by searching globally with different core compositions. 
The obvious next step is to extend this method to treat the cen- 
tral oxygen abundance as a free parameter. This has the exciting 
potential to place meaningful constraints on the 12 C(a,7) 16 
nuclear reaction cross-section. 
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Table 1 

Convergence of the method on synthetic data. 



Iteration 


?eff 


M*/M 


log(M He /M*) 


X 


q 


1" 


23,600 


0.600 


-5.76 


0.52 


0.55 


\ b 


22,200 


0.660 


-2.79 


0.99 


0.55 


2 a 


22,200 


0.660 


-2.79 


0.88 


0.51 


2 b 


22,600 


0.650 


-2.74 


0.85 


0.51 


3 a 


22,600 


0.650 


-2.74 


0.80 


0.50 


3 b 


22,600 


0.650 


-2.74 


0.80 


0.50 



" Value of M^/Mq fixed during this iteration. 
b Value of q fixed during this iteration. 



We do not presently include oxygen in the envelopes, so the 
mass fraction of oxygen must drop to zero by this point. 

We calculated the magnitude of deviations from the mean 
period spacing for models using our profiles compared to those 
due to smooth profiles from recent theoretical calculations by 
Salaris et al. (1997). The smooth theoretical profiles caused 
significantly larger deviations, so we conclude that the abrupt 
changes in the oxygen abundance resulting from our parame- 
terization do not have an unusually large effect on the period 
spacing (see Appendix for details). Although the actual chemi- 
cal profiles will almost certainly differ from the profiles result- 
ing from our simple parameterization, we should still be able 
to probe the gross features of the interior structure by matching 
one or another linear region of the presumably more compli- 
cated physical profiles. 

We used the same ranges and resolution for M„, T e ff, and Mee 
as in MNW, so the search space for this 5-parameter problem is 
10,000 times larger than for the 3-parameter case. Initially we 
tried to vary all 5 parameters simultaneously, but this proved 
to be impractical because of the parameter-degeneracy between 
M* and q. The pulsation properties of the models depend on the 
radial location in the chemical profile where the oxygen mass 
fraction begins to change. If the genetic algorithm (GA) finds 
a combination of and q that yields a reasonably good fit to 
the data, most changes to either one of them by itself will not 
improve the fit. As a consequence, simultaneous changes to 
both parameters are required to find a better fit, and since this is 
not very probable the GA must run for a very long time. Tests 
on synthetic data for the full 5-parameter problem yielded only 
a 10 percent probability of finding the input model even when 
we ran for 2000 generations — ten times longer than for the 3- 
parameter case. By contrast, when we used a fixed value of 
q and repeated the test with only 4 free parameters, the GA 
found the input model in only 400 generations for 8 out of 10 
runs. Even better, by fixing the mass and allowing q to vary, 
it took only 250 generations to find the input model in 7 out 
of 10 runs. This suggests that it might be more efficient to al- 
ternate between these two subsets of 4 parameters, fixing the 
fifth parameter each time to its best-fit value from the previous 
iteration, until the results of both fits are identical. 

Since we do not know a priori the precise mass of the white 
dwarf, we need to ensure that this iterative 4-parameter ap- 
proach will work even when the mass is initially fixed at an 
incorrect value. To test this, we calculated the pulsation pe- 
riods of the best-fit 0.65 M Q model for GD 358 from MNW 
[r eff = 22,600 K, log(M He /M*) = -2.74, 20:80 C/O "shallow"] 
and then iteratively applied the two 4-parameter fitting routines, 
starting with the mass fixed at 0.60 M Q — a difference compa- 



rable to the discrepancy between the mass found by MNW and 
the value found by Bradley & Winget (1994). The series of fits 
leading to the input model are shown in Table 1 . This method 
required only 3 iterations, and for each iteration we performed 
10 runs with different random initialization to yield a probabil- 
ity of finding the best-fit much greater than 99.9 percent. In 
the end, the method required a total of 2.5 x 10 6 model evalua- 
tions (128 trials per generation, 10 runs of 650 generations per 
iteration). This is about 200 times more efficient than calculat- 
ing the full grid in each iteration, and about 4,000 times more 
efficient than a grid of the entire 5 -dimensional space. 

Next, we applied this method to the observed pulsation pe- 
riods of GD 358. We initially fixed the mass at 0.61 M Q , the 
value inferred from the original asteroseismological study by 
Bradley & Winget (1994). The solution converged after four 
iterations, and the best-fit values of the five parameters were: 

T eB = 22, 600 K 

M*/M Q = 0.650 

log(M He /M*) = -2.74 

X = 0.84 

q = 0.49 

Note that the values of M„, T e ff and Mee are identical to those 
found in MNW. The best-fit mass and temperature still dif- 
fer significantly from the values inferred from spectroscopy by 
Beauchamp et al. (1999). However, the luminosity of our best- 
fit model is consistent with the luminosity derived from the 
measured parallax of GD 358 (Harrington et al. 1985). 

To alleviate any doubt that the GA had found the best com- 
bination of Xq and q, and to obtain a more accurate estimate 
of the uncertainties on these parameters, we calculated a grid 
of 10,000 models with the mass, temperature, and helium layer 
mass fixed at their best-fit values. A contour plot of this grid 
near the solution found by the GA is shown in Figure 1 . 

3. REVERSE APPROACH 

The results of forward modeling with one adjustable point 
in the chemical profile make it clear that information about the 
internal structure is contained in the data, and we just need to 
learn how to extract it. If we want to test more complicated ad- 
justable profiles, forward modeling quickly becomes too com- 
putationally expensive as the dimensionality of the search space 
increases. We need to devise a more efficient approach to ex- 
plore the myriad possibilities. 

The natural frequency that dominates the determination of 
pulsation periods in white dwarf models is the Brunt-Vaisala 
(B V) frequency, which we calculate using the Modified Ledoux 
treatment described in Tassoul Fontaine & Winget (1990). To 
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FIG. 1. — A contour plot of the central oxygen mass fraction (Xo) versus 
the fractional mass location of the change in the oxygen gradient (q), with M* , 
r e ff, and Mhc/M* fixed at their best-fit values. The model with the absolute 
minimum residuals (identical to the best-fit found by the GA) is marked with 
an X. The contours are drawn at 1, 3, 10, 25 and 40 times the observational 
noise. 

get a sense of how the pulsation periods depend on the BV 
frequency in various regions of the model interior, we added 
a smooth perturbation to the best-fit model of GD 358 from 
MNW, moving it one shell at a time from the center to the sur- 
face. The perturbation artificially decreased the BV frequency 
across seven shells, with a maximum amplitude of 10 percent. 
We monitored the effect on each of the pulsation periods as the 
perturbation moved outward through the interior of the model. 
The results of this experiment for the pulsation periods corre- 
sponding to those observed in GD 358 are shown in Figure 2. 

Essentially, this experiment demonstrates that the pulsation 
periods are sensitive to the conditions all the way down to the 
inner few percent of the model. Since the observational uncer- 
tainties on each period are typically only a few hundredths of a 
second, even small changes to the BV frequency in the model 
interior are significant. 

The root-mean-square (rms) residuals between the observed 
pulsation periods in GD 358 and those calculated for the best-fit 
from forward modeling are still much larger than the observa- 
tional noise. This suggests that either we have left something 
important out of our model, or we have neglected to optimize 
one or more of the parameters that could in principle yield a 
closer match to the observations. To investigate the latter possi- 
bility, we introduced ad hoc perturbations to the BV frequency 
of the best-fit model to see if the match could be improved. Ini- 
tially, we concentrated on the region of the Brunt-Vaisala curve 
that corresponds to the internal chemical profile. 

If we look at the BV frequency for models with the same 
mass, temperature, helium layer mass, and central oxygen mass 
fraction but different internal chemical profiles (see Figure 3) it 
becomes clear that the differences are localized. In general, we 
find that changes in the composition gradient cause shifts in the 
BV frequency. Moving across an interface where the gradient 
becomes steeper, the BV frequency shifts higher; at an interface 
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FIG. 2. — For the best-fit model of GD 358 from Metcalfe Nather & Winget 
(2000) this plot shows the change in pulsation period for 1=1 modes of various 
radial overtone number (k) which result from a smooth artificial 10% decrease 
in the Brunt-Vaisala frequency as a function of the natural log of the ratio of 
the distance from the center of the model to the local pressure (top axis) and 
the fractional mass coordinate -log(l -m/M„) (bottom axis). The center of 
the model is to the left, and the surface is to the right. Also indicated is the 
mass fraction expressed as a percentage for several values closer to the center. 




FIG. 3. — The Brunt-Vaisala frequency as a function of the radial coordi- 
nate ln(r/p) for several models with the same mass, temperature, helium layer 
mass, and central oxygen mass fraction but different internal chemical profiles 
(a) from the center of the model at left to the surface at right, and (b) only in 
the range of ln(r/p) indicated by the dashed box in the upper panel. The three 
curves correspond to a profile with q equal to 0.2 (dotted), 0.49 (solid), and 0.8 
(dashed). 
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where the gradient becomes more shallow, the BV frequency 
shifts lower. The greater the change in the gradient, the larger 
the shift in the BV frequency. 

3.1. Proof of Concept 

We began by generating a model with the same mass, tem- 
perature, helium layer mass, and internal composition as the 
best-fit from MNW, but using a uniform internal chemical pro- 
file with constant 20:80 C/O out to the 95 percent mass point. 
We identified a sequence of 100 shells in this model spanning a 
range of fractional mass from 0.20 to 0.97 and perturbed the B V 
frequency to try to produce a better match to the observations. 
We parameterized the perturbation as a multiplicative factor ap- 
plied to the BV frequency over a range of shells, described by 
four parameters: (1) the innermost shell to perturb, (2) the mag- 
nitude of the perturbation at the innermost shell, (3) the number 
of shells in the perturbation range, and (4) the magnitude of 
the perturbation at the outermost shell. These 4 parameters are 
sufficient to describe a profile with two abrupt changes in the 
composition gradient. 

The innermost shell was allowed to be any of the 100 shells, 
and the number of shells in the perturbation range was allowed 
to be between and 100. If the chosen range would introduce 
perturbations outside of the sequence of 100 shells, the outer- 
most of these shells was the last to be perturbed. The mag- 
nitude of the perturbation was allowed to be a multiplicative 
factor between 1.0 and 3.0 at both the innermost and outermost 
perturbed shells, and was interpolated linearly in the shells be- 
tween them. 

By using a GA to optimize the perturbation, we can try many 
random possibilities and eventually find the global best-fit with 
far fewer model evaluations than a full grid of the parameter 
space. We demonstrated this by introducing a particular per- 
turbation to the model and determining the pulsation periods. 
Using the same unperturbed model, we then attempted to find 
the parameters of the perturbation by matching the periods us- 
ing the GA. In 9 out of 10 runs (500 generations of 64 trials), 
the GA found a first-order solution within two grid points of 
the input perturbation. Thus, the GA combined with a small 
(625 point) grid search yields a 90 percent probability of suc- 
cess for an individual run. By repeating the process several 
times, the probability of finding the correct answer quickly ex- 
ceeds 99.9 percent, even while the number of model evaluations 
required remains hundreds of times lower than a full grid of the 
parameter- space. 

3.2. Application to GD 358 

Having demonstrated that the method works on calculated 
model periods, we applied it to the observed pulsation periods 
of GD 358. We began with a model similar to the new best-fit 
determined from the forward modeling described in §??, but 
again using a uniform internal chemical profile (constant 16:84 
C/O out to 0.95 m/Af*). After the GA had found the best-fit per- 
turbation for GD 358, we reverse-engineered the corresponding 
chemical profile. 

To accomplish this, we first looked in the unperturbed model 
for the fractional mass corresponding to the innermost and out- 
ermost shells in the perturbation range. We fixed the oxygen 
abundance to that of the unperturbed model from the center out 
to the fractional mass of the innermost perturbed shell. The 
size of the shift in the BV frequency is determined by how 
much the composition gradient changes at this point, so we ad- 



justed the oxygen abundance at the fractional mass of the outer- 
most perturbed shell until the change in the gradient produced 
the required shift. Finally, we fixed the oxygen abundance to 
that value from the outermost perturbed shell out to a fractional 
mass of 0.95, where it abruptly goes to zero. 

After we found the C/O profile of the best-fit perturbation in 
this way, we fixed this reverse-engineered profile in the mod- 
els and performed a new fit from forward modeling with the 
GA to re-optimize the mass, temperature, helium layer mass, 
and central oxygen mass fraction. The BV curve of the final 
model differs slightly, of course, from that of the original uni- 
form composition model with the perturbation added. But the 
approximate internal structure is preserved, and leads to a bet- 
ter match to the observed pulsation periods than we could have 
otherwise found. 



4. RESULTS 

The calculated periods and period spacings (AP = Pf,+i -Pk) 
for the best-fit models from the new forward modeling and from 
the reverse approach are shown in the bottom two panels of Fig- 
ure 4 along with the data for GD 358. The best-fit models of 
Bradley & Winget (1994) and MNW are shown in the top two 
panels for comparison. The data in Figure 4 for the observa- 
tions and our best-fit models are given in Table 1 . Some of the 
improvement evident in the panels of Figure 4 is certainly due 
to the fact that we have increased the number of free parame- 
ters. To evaluate whether or not the new models represent a 
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FIG. 4. — The periods and period spacings observed in GD 358 (solid points) 
with the theoretical best-fit models (open points) from (a) Bradley & Winget 
(1994), (b) Metcalfe Nather & Winget (2000), (c) the new forward modeling 
presented in §??, and (d) the reverse approach presented in §??. Uncertainties 
on the observations are smaller than the size of the points in this figure. 
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Table 2 

Periods and Period Spacings for GD 358 and the Best-fit Models. 





Observed 


3-parfit 


5-par 


fit 


7-par fit 


k 


P 


AP 


P 


AP 


P 


AP 


P 


AP 


8 . . . 


423.27 


40.96 


422.31 


42.26 


422.36 


40.92 


422.75 


39.69 


9... 


464.23 


37.36 


464.57 


36.77 


463.28 


38.01 


462.43 


37.46 


10... 


501.59 


40.16 


501.35 


35.88 


501.29 


37.50 


499.90 


39.70 


1 1 

ii... 


541.75 


35.01 


537.23 


39.04 


7Q 
J JO. 1 y 


^f, 77 


539.60 


36.65 


12... 


576.76 


41.52 


576.27 


42.79 


575.56 


43.33 


576.25 


42.10 


13... 


618.28 


40.07 


619.06 


39.79 


618.89 


39.85 


618.36 


40.25 


14... 


658.35 


42.29 


658.85 


42.97 


658.74 


42.36 


658.61 


42.90 


15... 


700.64 


33.66 


701.82 


32.76 


701.10 


34.33 


701.51 


33.44 


16... 


734.30 


36.37 


734.58 


36.92 


735.42 


36.99 


734.95 


36.59 


17... 


770.67 


40.03 


771.50 


39.30 


772.41 


37.34 


771.54 


39.60 


18... 


810.7 


44.1 


810.80 


44.34 


809.75 


44.63 


811.14 


44.15 



significant improvement to the fit, we use the Bayes Infor- 
mation Criterion (BIC), following Montgomery Metcalfe & 
Winget (2001). 

The fit in MNW used n p = 3 completely free parameters, 
though they also sampled several combinations of two addi- 
tional parameters. This amounts to a partial optimization in 5 
dimensions. To make a fair prediction of the residuals required 
to be considered significant we use their best-fit carbon core 
model, which represents the best truly 3-parameter fit. This 
model had rms residuals of a(P) = 2.30 seconds for the peri- 
ods and cr(AP) = 2.65 seconds for the period spacings. For 
N = 1 1 data points, the BIC leads us to expect the residuals of 
a n p = 5 fit to decrease to a{P) = 1.84 and a(AP) = 2.13 just 
from the addition of the extra parameters. In fact, the fit from 
the new forward modeling presented in §?? has a(P) = 1 .28 and 
<j( AP) = 1.42, so we conclude that the improvement is statisti- 
cally significant. 

The results of the reverse approach presented in §?? are 
harder to evaluate because we are perturbing the BV frequency 
directly, rather than through a specific parameter. We consider 
each additional point in the internal chemical profile where the 
composition gradient changes to be a free parameter. Under this 
definition, the perturbed models are equivalent to a 7-parameter 
fit since there are three such points in the profiles, compared to 
only one for the 5-parameter case. If we again use the BIC, 
we expect the residuals to decrease from their n p = 5 values 
to a(P) = 1.03 and a(AP) =1.14 seconds. After re-optimizing 
the other four parameters using the profile inferred from the re- 
verse approach, the residuals actually decreased to a(P) =1.11 
and <r(AP) = 0.71 seconds. The decrease in the period residuals 
is not significant, but the period spacings are improved consid- 
erably. This is evident in the bottom panel of Figure 4. 

Adopting the central oxygen mass fraction from the best-fit 
forward modeling, Xo = 0.84 ± 3 we can place preliminary con- 
straints on the 12 C(a,7) 16 cross-section. Salaris et al. (1997) 
made detailed evolutionary calculations for main-sequence stel- 
lar models with masses between 1 and 7 M Q to provide internal 
chemical profiles for the resulting white dwarfs. For the bulk of 
the calculations they adopted the rate of Caughlan et al. (1985) 
for the 12 C(a,7) 16 reaction (S 30 o = 240 keV barns), but they 
also computed an evolutionary sequence using the lower cross- 
section inferred by Woosley Timmes & Weaver (1993) from 
solar abundances (S300 = 170 ±50 keV barns). The chemical 
profiles from both rates had the same general shape, but the 
oxygen abundances were uniformly smaller for the lower rate. 



In both cases the C/O ratio was constant out to the 50 percent 
mass point, a region easily probed by white dwarf pulsations. 

The central oxygen mass fraction is lower in higher mass 
white dwarf models. The rate of the triple-a reaction (a three- 
body process) increases faster at higher densities than does the 
12 C(a,7) 16 reaction. As a consequence, more helium is used 
up in the production of carbon, and relatively less is available to 
produce oxygen in higher mass models. Interpolating between 
the models of Salaris et al. (1997) which used the higher value 
of the cross-section, we expect a central oxygen mass fraction 
for a M = 0.65 M Q model of X^ sh = 0.75. Using additional 
calculations for the low rate (M. Salaris 2001, private commu- 
nication), the expected value is X ow = 0.62. Extrapolating to 
the value inferred from our new forward modeling, we estimate 
that the astrophysical S-factor at 300 keV for the 12 C(a,7) 16 
cross-section is in the range S300 = 290 ± 15 keV barns (internal 
uncertainty only). 

The internal chemical profiles corresponding to the best-fit 
models from §?? and §?? are shown in Figure 5 with the the- 
oretical profile for a 0.61 M Q model from Salaris et al. (1997), 
scaled to a central oxygen mass fraction of 0.80. The profile 
from the best-fit forward modeling matches the location and 
slope of the initial shallow decrease in the theoretical profile. 
The reverse approach also finds significant structure in this re- 
gion of the model, and is qualitatively similar to the Salaris et 
al. (1997) profile to the extent that our parameterization allows. 

5. DISCUSSION & FUTURE WORK 

The extension of the genetic-algorithm-based approach to 
optimize the internal composition and structure of our best-fit 
white dwarf models to GD 358 has yielded some exciting re- 
sults. The values of the 3 parameters considered in the initial 
study (M*, r e ff,A^He) are unchanged in the full 5-parameter fit, 
so we feel confident that they are the most important for match- 
ing the gross period structure. The significant improvement to 
the fit made possible by including Xo and q as free parameters 
confirms that the observed pulsations really do contain infor- 
mation about the hidden interiors of these stars. 

The efficiency of the GA relative to a grid search is much 
higher for the larger parameter-space, and the ability of the 
method to find the global best-fit is undiminished. The applica- 
tion of this method to data on additional pulsating white dwarfs 
will help us to assess the significance of our results, so we can 
begin to understand the statistical properties of these ubiquitous 
and relatively simple stellar objects. It may also provide us with 
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ln(r/p) 

FIG. 5. — The internal oxygen profiles (top) and the corresponding region 
of the Brunt- Vaisala' curves (bottom) for the best-fit forward model from §?? 
(dotted), the result of the best-fit reverse approach from §?? (dashed), and the 
scaled theoretical calculations of Salaris et al. (1997) for comparison (solid). 

new insights into the processes involved in white dwarf evolu- 
tion. 

Our reverse approach to model-fitting was originally moti- 
vated by the fact that the variance of our best-fit model from the 
forward method was still far larger than the observational un- 
certainties. Although this initial application has demonstrated 
the clear potential of this approach to yield better fits to the 
data, the improvement to the residuals was only marginally sig- 
nificant. We will continue to develop this approach, but we 
must simultaneously work to improve the input physics of our 
models. In particular, since the internal profile from this ini- 
tial application is qualitatively similar to the theoretical profiles 
(given what is possible with this parameterization), we should 
really include oxygen in our envelopes and eventually calculate 
fully self-consistent models out to the surface. 

Turning the central oxygen abundance into a more precise 
constraint on the 12 C(a,7) 16 nuclear reaction cross-section 
will require additional detailed simulations like those of Salaris 



et al. (1997). By determining the range of values for the cross- 
section that produce a central oxygen abundance within the 
measurement uncertainties of Xo, we should be able to sur- 
pass the precision of the extrapolation from laboratory data by 
nearly an order of magnitude. To assess the systematic uncer- 
tainties, we will also need to determine the relative importance 
of model-dependent details such as convection. 

We would like to thank Ed Nather and Mike Montgomery for 
helpful discussions, and M. Salaris for providing us with data 
files of the theoretical white dwarf internal chemical profiles. 
We are grateful to the High Altitude Observatory Visiting Sci- 
entist Program for fostering this collaboration in a very produc- 
tive environment for two months during the summer of 2000. 
This work was supported by grant NAG5-9321 from the Ap- 
plied Information Systems Research Program of the National 
Aeronautics & Space Administration, and in part by grant AST- 
9876730 from the National Science Foundation. 



: 1 i 1 1 1 i 1 1 1 i 1 1 1 r 




-34 -32 -30 

ln(r/p) 





50 


_ 1 




,1111,1 








\ 


/\ / 










\ 






x /' / 


1 






\ y\ \ t 




X 7 / 

\'7 7 




40 




vA /I \ \ ■■ 


VS^.\ ,'k 


a. 


30 


" 1 


>v6~o..\ <[' 

\ i v° \f 
v \ / 

, i , , 


" X_ ">/ 


v° / / 
\ / 7 

, i i , , , - 






400 


500 


600 700 
Period (sec) 


800 



FIG. 6. — The Brunt- Vaisala frequency in the region of the models that 
correspond to the internal chemical profile (top panel) and the period spacing 
diagrams (bottom panel) for our best-fit forward model from §2 (dotted), a uni- 
form 20:80 C/O model (dashed), and a scaled theoretical profile from Salaris 
et al. (1997) for comparison (solid). 



APPENDIX 

At first glance, it may seem that our parameterization of the internal chemical profile might introduce large deviations from the 
mean period spacing due to the abrupt change in the composition gradient. To demonstrate that the models are not unusually sensitive 
to our profiles when compared to smooth profiles, we calculated the pulsation periods of models with the best-fit parameters from 
Metcalfe Nather & Winget (2000) using several internal chemical profiles. The results are shown in Figure 6. In the top panel, 
we show the Brunt-Vaisala frequency in the region of the models that correspond to the internal chemical profile. In the bottom 
panel, we show the period spacing diagrams for each of these models. We find that deviations from the mean period spacing for a 
model using an internal chemical profile from our parameterization have a magnitude comparable to those caused by a profile with a 
uniform C/O mixture out the the 95 percent mass point. The smooth profile from Salaris et al. (1997) causes the largest deviations of 
the three. 
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